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The tensor-virial method is appUed for a study of osciUation modes of uniformly rotating Bose- 
Einstein condensed gases, whose rigid body rotation is supported by an vortex array. The second 
order virial equations are derived in the hydrodynamic regime for an arbitrary external harmonic 
trapping potential assuming that the condensate is a superfluid at zero temperature. The axisym- 
metric equilibrium shape of the condensate is determined as a function of the deformation of the 
' trap; its domain of stability is bounded by the constraint f2 < 1 on the rotation rate (measured 

' in units of the trap frequency ujq.) The oscillations of the axisymmetric condensate are stable 

. with respect to the transverse-shear and toroidal modes of oscillations, corresponding to the / = 2, 

^ ' \m\ — 1,2 surface deformations. The eigenfrequencies of the modes are real and represent undamped 

, oscillations. The condensate is also stable against quasi-radial pulsation modes (Z = 2, m = 0), and 

■ its oscillations are undamped, if the superflow is assumed incompressible. In the compressible 

case we find that for a polytropic equation of state, the quasi-radial oscillations are unstable when 
7(3 - n^) < 1 - 30.'^, and are stable otherwise. Thus, a dilute Bose gas, whose equation of state is 
polytropic with 7 = 2 to leading order in the diluteness parameter, is stable irrespective of the ro- 
' tation rate. In non-axisymmetric traps, the equilibrium constrains the (dimensionless) deformation 

^ , in the plane orthogonal to the rotation to the domain A2 > fl^ with Q, < 1. The second harmonic 

' oscillation modes in non-axisymmetric traps separate into two classes which have even or odd parity 

0^ I with respect to the direction of the rotation axis. Numerical solutions show that these modes are 

CO . stable in the parameter domain where equilibrium figures exist. 
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I. INTRODUCTION 



' After the experimental realization of Bose-Einstein condensation in vapors of alkali atoms , the understanding of 
their behavior under rotation became the main focus of theoretical and experimental work. Experimental observation 
of the vortex states in a two-component system Q , and the realization of the vortex states in a stirred one-component 
' O Bose-Einstein condensate |^ confirmed the expectations that sufficiently large condensates exhibit bulk superfluid 
gjperties by supporting their rotation above the critical angular frequency by the Fcynman-Onsager vortex state 

The theoretical approaches to the Bose-Einstein condensed gases commonly distinguish between the two regimes 
of the strong and the weak interparticle interactions, which correspond to the conditions Na/d 3> I and Na/d <C 1, 
respectively, where N is the number of particles in the gas, a is the scattering length, and d — y/Ti/mujo is the oscillator 
^ ' length defined in terms of the oscillator frequency and the boson mass m (for a review see e.g. []8|,|lT|.) In the weak 
coupling regime, where the scale of the variation of the condensate wave-function is of the same order as the size of 
the cloud, quantum effects play an essential role and the system behaves much like an atomic nucleus (see e. g. Ref. 
[ p2[ ) . In the strong coupling regime and for sufficiently large condensate, the condensate is well described in terms of 
hydrodynamics. Under these conditions the time-dependent Gross-Pitaevskii (G-P) equation for the (complex) 
condensate wave function reduces to two (real) hydrodynamic equations for the density and velocity of the condensate 
at zero temperature. One consequence of the assumption that the system is sufficiently large is that the quantum 
ressure term can be ignored compared with the interaction terms. Another is that the Thomas-Fermi approximation 
is valid and the equilibrium kinetic energy density is negligible compared to the potential energy density. This 
approximation also ensures that the thermodynamic quantities like the pressure or the chemical potential of the gas 
are well defined at any local point of the gas cloud. Finally, the coherence length of the condensate ^ = l/\/8nna, with 
n being the number density of condensate particles, becomes sufficiently small in this limit to allow for a description 
of a vortex core as a singularity. 

The superfluid hydrodynamic approach has been used in the studies of the collective excitations of the condensate 
clouds [||jl5| the scissors mode oscillations , the moment of inertia of a rotating condensate , mainly in the 
irrotational regime; see however Ref. |Q. The critical angular velocity of condensate, fid, at which the creation of a 
vortex is energetically favorable, scales, in units of the trap frequency, as ilci oc (d^/i?^) In (i?/^), where R is the size 
of the cloud; the numerical coefficient depends on the details of the geometry of the trap 1^. For a sufficiently large 
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condensate the ground state would correspond to coarse-grained rigid-body rotation supported by a vortex lattice. In 
this case, the irrotational constraint on the average superfluid velocity V X v = is replaced by the rigid-body rotation 
condition Vxv = 2J7. The purpose of this paper is the study of the oscillations of the trapped Bose condensates 
under rigid-body rotation. We shall not discuss the modes related to the vortex structure itself (e.g. Tkachenko 
modes). We shall confine ourselves to the zero-temperature limit, in which case the effects of mutual friction due to 
the interaction of the thermal excitations with the vortex state, and dissipation due to the viscosity of the thermal 
component both can be neglected. 

Below, we shall apply the tensor virial method originally developed for the study of equilibrium and stability of 
rotating liquid masses bound by self-gravitation The tensor virial method transforms the local hydrodynamical 
equations into global virial equations that contain the full information on the structure and stability of a system 
as a whole. The method is especially useful for studying perturbations of incompressible uniform ellipsoids from 
equilibrium, in which case the each perturbed virial equation yields (in the absence of viscous dissipation) a different 
set of normal modes. The extension to compressible flows is straightforward in the case of adiabatic perturbations. 
Moreover, when the equilibrium is sustained by an external confining potential, fluid perturbations do not back-react 
(change) the confining potential itself, and the tensor virial method can be extended to non-uniform compressible 
flows. The description is particularly useful for gases with polytropic equations of state p oc p'^ , where p is the 
pressure, p is the density, and 7 is the adiabatic index. In that case, we shall see that all of the modes of the gas 
can be found from the tensor virial equations. (We illustrate this only for the second harmonic modes of oscillation 
in this paper, but the extension to higher order modes is straightforward, and will be discussed elsewhere.) For an 
interacting Bose gas the pressure is a non- analytic function of the diluteness parameter na^. However the equation 
of state to the leading order in the small parameter na^ <C 1 can be written in a polytropic form: p = Kp^, with 
Since the zero temperature equations of motion of a trapped rotating condensate turn out formally 
identical to the corresponding equations of motion of a ordinary (non-superfluid) liquid, our results might be of 
significance in a broader context. On the other hand, present results can serve as a starting point for an extension to 
finite temperatures where, in addition to the superfluid, the fluid of normal quasiparticle excitations plays a role. 

The paper is organized as follows. In § 2 we derive the first and second order virial equations for trapped Bose- 
condensed gases in the hydrodynamic regime. The equilibrium shape of rotating condensates and the second harmonic 
modes of oscillations in axisymmetric traps are discussed in § 3. § 4 discussed the second harmonic modes of 
oscillations in non-axisymmetric traps. Our results are summarized in § 5. The Appendix gives a brief derivation of 
the hydrodynamic equations from the G-P equation. 



II. VIRIAL EQUATIONS FOR A TRAPPED CONDENSATE 

The Euler equation for a condensate in a harmonic, arbitrary deformed trap, can be written as 

(d d\ dp pdK , p 9|»xxp , ^ 

' U ^ "9^ ~ 2 + + 2pe.„..,a., (1) 

where the Latin subscripts denote coordinate directions; p, p, and u are the density, pressure, and velocity of conden- 
sate; (summation over repeated indexes is assumed). The external harmonic trapping potential is 

3 

0tr(x)=^o'E^'^?' (2) 

2=1 

where Aj are the dimensionless deformation parameters, Wo is the frequency of the harmonic oscillator in the trap. 
Eq. (^ is written in a frame rotating with angular velocity relative to some inertial coordinate reference system. 
Its derivation from the G-P equation for the condensate wave-function is given in the Appendix. Our fundamental 
assumption is that the condensate undergoes a rigid-body rotation, which is supported by an array of vortices. In this 
case, as explained in the Appendix, the equations of motion are formally identical to those of a ordinary compressible 
(non-superfluid) liquid. 

Starting from Eq. (|^) we can construct a set of virial equations of various order. Suppose that the condensate 
occupies a volume V . Taking the zeroth moment of Eq. (|l|) amounts to integrating over V\ doing so, we obtain 

d^xpui^ = 2eii,n^.m J d^xpui + {il'^Sij - Jl,^ ) J d^xpxj 
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The first order virial equation describes the center of mass motion of the condensate, and is trivial as these motions 
correspond to a uniform translation of the system as a whole. 

Taking the first moment of Eq. results in the second order virial equation 

+2T,j + 5,jIV, (4) 

where 

lij = I (fixpxiXj, 11=1 (fixp, Tij = — I (fixpuiUj. (5) 
Jv Jv 2 Jy 

Consider the variation of the second order virial equation under the influence of perturbations. The Eulerian 
variations of various terms in Eq. are straightforward |20| and we find 

i d?xP'^i^] = (^^ ^ ^lAi)Vij - Q.iQ.kVk] + 2e,;;mfi„i5 / tfxpuiXj + 2STij + SijSU, (6) 
Jv Jv 

where for a Lagrangian displacement ^ 

Vij=6Iij= / (fxp{^iXj + (,jXi). (7) 
Jv 

The tensors Vij are manifestly symmetric in their indexes, as is evident from their definition. It is useful to define 
also their non-symmetric parts as 



V^■J = / d^xp£,^Xj, (8) 

V 



such that Vij = Vi-j + Vj-i. When there are no condensate motions in the unperturbed state in the rotating frame 
(u = 0), the second order virial equation becomes 



The tensor virial Eq. can be used for the study of small amplitude oscillations of a trapped condensate. For 
time-dependent Lagrangian displacements of the form 

i{x,,t)^i{x,)^-' (10) 

Eq. becomes 

\^V,,j - 2^eaA^m,, = - i^lAi)V^j - ^^nkVkJ + %(5n. (11) 
Eq. (O) contains all the second harmonic modes of the rotating condensate in a arbitrary harmonic trap. 



III. EQUILIBRIUM AND OSCILLATIONS IN AXISYMMETRIC TRAPS 



A. Equilibrium shape 

Next, we specialize Eq. (^) to the case of axisymmetric traps, Ai = A2 ^ A3 and assume, without loss of generality, 
Ai = 1. The assumption of the axisymmetry is, in practice, an approximation as a perfectly axisymmetric trap will 
not exert a torque on the condensate. Our assumption is that the deviations from the axisymmetry, which drive the 
rotation of the condensate, are small compared with the deformation of the trap \ Ai — A2I/A1 <C 1. In principle, once 
the condensate is brought to rotation the axisymmetry can be restored. 

Consider a equilibrium state in which the condensate rotates uniformly with the rotation vector along the X3 axis 
of Cartesian system of coordinates. Then Eq. (Q) reduces to 

{hj - Sahj) - uolA.Uj = -S,jU. (12) 
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In axisymmetric traps with 03 denoting the eUipsoidal semi-axis along the rotation, the two remaining semi-axis of 
the elhpsoid, ai and a2, are equal. The diagonal component of Eq. (|l2| ) provides the relation between the rotation 
frequency, the deformation parameters along the direction (both are fixed in an experiment) and the ratio of the 
semi-axis of the resulting spheroidal figure: 



4^:^(1-2,1. (13) 



Depending whether 03 < oi or 03 > ai the condensate shape is either oblate or prolate. Note that the rotation 
frequency is bounded from above, Q/ujq < 1, and for A3 > 1 the equilibrium figures are always oblate. A striking 
consequence of the Eq. (|l3|), when A3 = 1 — il^/cjQ, is an equilibrium figure which is a rotating sphere Although 
we had to assume small deviations from axisymmetry to make the condensate rotate, this assumption is a prerequisite 
for producing rotation, but not an intrinsic property of the rotating condensate. 

The equilibrium density profile of the condensate, for a polytropic equation of state, can be obtained by a direct 
integration of the unperturbed limit of Eq. (|l|): 

p{^)=Po\l--^^^[{l-^')xi+{A2-n^)xl + A3xl]\ e[...], (14) 

where po is the central density of the cloud, the dots in the argument of the 0-function (which insures that the 
condensate density is positive) stand for the expression in the curly brackets. Here we have retained the parameter 
A2 Ai to incorporate the case of non-axisymmetric traps which is treated in the next section. The effect of the 
rotation is the stretching of the profile of the condensate in the equatorial plane due to the centrifugal potential. 
According to Eq. ([l^), the density is constant on concentric ellipsoids in the unperturbed rotating background. Note 
that the normalization of the wave-function of the condensate implies 

d^xp{x)^M, (15) 

V 

where M is the total mass in the cloud. 



B. Second harmonic modes of oscillation 



Next consider perturbations from the equilibrium state of uniform rotation, with the spin vector along the X3 axis. 
Surface deformations related to various modes can be classified by corresponding terms of the expansion in surface 
harmonics labeled by indexes l,m. Second order harmonic deformations correspond to I — 2 with five distinct values 
of TO, —2 < TO < 2. The 18 equations represented by Eq. ( pi] ) separate into two independent subsets which are odd 
and even with respect to the index 3. The corresponding oscillation modes can be treated separately. 



1. Transverse-shear modes 



These modes correspond to the surface deformations with |to| = 1 and represent relative shearing of the northern 
and southern hemispheres of the spheroid. The components of Eq. (Kw, which are odd in index 3, are: 



X''V3-i = -A3V31, (16) 

AV3;2 = -A3V32, (17) 

\^Vi,3 ~ 21]AV^2;3 = -^13 + f^Vi3, (18) 

AV2;3 + 2nxvi,3 = -V23 + n^V23. (19) 



We sum the Eqs. ([I6D, ( |18[ ) and (|17|), (|19D, respectively, and use the symmetry properties of Vij combined with Eqs. 
(^, (^. This results in the relations 

A (A^ + A3 + 1 - Vi3 - 2fl(A2 + A3)V23 = 0, (20) 
A(A2 +A3 + I- n^)V23 - 2f^(A2 + ^3)^23 = 0. (21) 
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The characteristic equation can be factorized by substituting \ = ia and we find 

- 2^20-2 _ (1 + ^3 _ fj2y ^ 2rjA3 = 0. 



The roots are give by 
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2n 1 



f2.3 



(,s+ + ,s_) ± 



(22) 



(23) 



where 



1 - 2A, 
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1 + A3 



1/2 



(24) 



Three complementary modes fohow from Eqs. (^3|)-(^^ via the replacement H, —fl. It can be verified that the 
condition 



n 



- 1-2A.- — 



n 2 



1 



- ( 1 + A3 + — 



n 3 



<o, 



(25) 



is satisfied for any A3 and < 1, therefore all three roots are real. The real frequencies of the transverse-shear 
modes are shown in Fig. 1 as a function of A3. Note that in the zero-temperature limit, to which present analysis is 
restricted, these modes are purely real, i. e. represent undamped oscillations. 



2. Toroidal modes 

These modes correspond to |m| = 2 and the motions in this case are confined to the planes parallel to the equatorial 
plane. The components of Eq. (pi]), which are even in index 3, are: 

X^V3.,3 = SI1-A3V33, (26) 

X'^Vi-i - 2nXV2;i = 611 - l)Vn, (27) 

AV2;2 + 2f7AFi;2 =611+ - 1)^22, (28) 

AVi;2-2r!AT/2;2-(f^'-l)T42, (29) 

AV2;i + 2nXVia = (^^^ - 1)^^21- (30) 



We add Eqs. (29) and (|30|), and subtract Eqs. ( pTI) and ( |28| ) to find the following coupled equations: 

[A^ + 2(1 - n^)] {Vn - V22) - mXVu = 0, (31) 
[A^ + 2(1 - n^)] V12 + nX{Vii - V22) = 0. (32) 

The characteristic equation for the toroidal modes is 

[A^ - 2{n^ - 1)] ' + ix'^n'^ = 0, (33) 

which is factorized by writing A = ia. Two solutions are then 



ai,2 = n± V2^^. (34) 

There are two complementary modes which are found by substituting —fl for Q. Since the rotation frequency is 
bounded (17 < 1) the toroidal modes are stable independent of the magnitude of the deformation in the equatorial 
plane. For the same reason these modes are undamped in the zero-temperature limit considered here. 
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3. Pulsation modes 



To find the pulsation modes, which correspond to m — 0, we first add the Eqs. (P7|)-(P^), and then use the Eq. 
( p6| ) to eliminate SIl in the result. In this manner we find that 

(AV2 + {Vn + V22) + 2nX{Vi,2 - V2a) - (A" + 2^3)^33 = 0. (35) 
Subtracting Eqs. ( |30| ) and ( p9|) (and discarding the A = root) one finds 

A(yi;2-V^2;l)-f^(^ll+V^22)=0. (36) 



Eqs. ( |35| ) and (36) can be further combined to a single equation: 

(A^ + 21^2 + 2) [Vii + V22) - 2(A2 + 2^3)^33 = 0. (37) 



It is instructive first to consider the case where the superfluid is incompressible. Then the Eq. (|3^) should be 
supplemented by the divergence free condition 

Yii±}:ji + % = 0. (38) 

Again, writing X ~ ia, we find for the square of frequency of the pulsation mode in the incompressible limit 

As CTg is always positive, these modes correspond to undamped stable oscillations. 

For compressible fluids we need the variation of the pressure tensor, which for adiabatic perturbations can be 
written as 

m = (7 - 1) Jd^x V,p, (40) 

where we assumed a polytropic equation of state p — Kp"' ; the polytropic index for a Bose gas to leading order in the 
parameter pa^ is equal 2. To evaluate the gradient of the pressure we turn to the Euler Eq. (|^) in the unperturbed 
limit, 

dh I dp _ ld(j)tr ^ I 9|»xxp ^^^^ 
dxi p dxi 2 dxi 2 dxi ' 

where the enthalpy h is defined by dh = dp/ p and is simply h — p^^^ / — 1) for p — Kp'' . Substituting the 
explicit expression for the trapping potential, we find that the gradient of the enthalpy is a linear function of the 
coordinates and, hence, the variation of the pressure tensor can be expressed in terms of virials Vu: 

= -t^[{y^^ + V22){1 - n") + ^3^33]. (42) 

Any of the equations which are even in the index 3 now can be used to close the system of equations for V33 and the 
virial combination Vn + V22- Substituting Eq. (42) for Sli in, e.g, Eq. ([2^), one finds 

[A2 + (7 + 1)A3]F33 + (7 - 1)(1 - ^""Wll + V22) = 0. (43) 

Eqs. (^^ and (^) completely determine the unknown virials; the characteristic equation for the pulsation modes, 
which is quadratic in A^, is 

A'^ + A^[7(A3 + 2 - 21^2) + A3 + An^] + 2yl3[7(3 - n'^) + 3f^2 - 1] = 0. (44) 

On substituting A = iti, the solution of the resulting quadratic equation becomes 

4 = i [7(^3 + 2 - 2^(2) + A3 + 4^*2] 

± ^ V[7(^3 + 2 - 2r!2) + A3 + 4f|2]2 _ 8A3[7(3 - Vf^) + - 1]. (45) 
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It is easy to see that there are only unstable modes if 7 < (1 — 3f7^)/(3 — fl"^) < 1/3; otherwise all modes are stable. 
In particular, all modes are stable for 7 = 2. Note, too, that there are twice as many solutions as were found in the 
incompressible case; Eq. ( ^9| ) can be recovered from Eq. ( p5| ) by taking the 7 — > 00 limit of cr^ . The origin of the 
additional root, cr^, may be traced to the dependence of Eq. ( [43|) , which only reduces to the incompressibility 
condition, Eq. (^8|), in the 7 ^ 00 limit if it is assumed that |A^| = |?| <C 7. 

The square of the pulsation modes as a function of the A3 is plotted in Fig. 3 for 7 = 2, the polytropic index of a 
dilute Bose gas. 



IV. EQUILIBRIUM AND OSCILLATIONS IN NON-AXISYMMETRIC TRAPS 



A. Equilibrium shape 



When the symmetry with respect to the rotation axis is broken, the equilibrium constraint Eq. (y_3|) needs to 
be supplemented by a relation fixing the semi-axis ratios in the plane orthogonal to the spin axis. The diagonal 
components of Eq. ( [l^ ) provide the triangle relations which determine the non-axisymmetric equilibrium figure 

al{n^ - LuiAi) = al{n^ - LjiA2) = -a^LuiAs. (46) 

These simultaneous constraints can be written in a equivalent form 



V ujf, J aj A3 \ uj^ 



Note that in addition to the upper bound on the rotation frequency set by the first relation (as in the case of the 
axisymmetric traps), the second relations places a lower bound on the deformation in the plane orthogonal to the 
spin axis: A2 > fP/ujQ. Given the experimentally controlled values of fl, A3 and A2, relations ( p7| ) determine, in a 
unique manner, the semi-axis ratios of the resulting figure. 



B. Second harmonic modes of oscillation 



The non-axisymmetric modes can be found from Eq. (y_l|) in a manner similar to the axisymmetric modes; however, 
now the degeneracy in indexes 1 and 2 should be relaxed. The oscillation modes separate into two non-combining 
groups which have even or odd parity with respect to the index 3. Below, we shall treat these modes separately. 



1. Odd modes 



Among the four components of Eq. (pT|), which are odd in index 3, three are identical to Eqs. (|T^), ([l7|), and ( |T^ ) 
under non-axisymmetric conditions; the component which is modifed reads 

AV2;3 + 217AFi;3 = + n^V23. (48) 

Summing Eqs. (|l|), (|l|), and (||), we arrive at 

A (A^ + A3 + 1 - 02) Vi3 - 2n{\^ + A3)V23 = 0, (49) 
A(A2 + A3 + A2- n^)V23 - 2f^(A2 + ^3)^23 = 0. (50) 

The sixth order characteristic equation derived from this algebraic system is 

A^ + [1 + ^2 + 2(^3 + n^)] \* + [A2 + A3+ A2A3 + Al-{1 + A2- 6A3)n^ + n^] + AAln' = 0. (51) 

The corresponding modes appear as conjugate pairs, i.e. there are only 3 distinct modes. The real parts of these 
modes (the modes are purely real) are show in Fig. 3 as a function of A3 for several values of O and fixed ratio 
r^VAs = 0.1. 
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2. Even modes 



The even parity components of the Eq. (11) which are given by the Eqs. (pq), (p7|), and ( p9[ ) remain unchanged 
when axisymmetry is relaxed; the remainder equations read 

x^V2;2 + 2nxvi,2 = sn + {n^ - A2)V2 



22, 



X'^V2-i + 2QXVi.i = {Q^ - A2)V2i. 



In the incompressible limit these equations should be supplemented by the condition 
of the second term) . In the compressible case the variations of the pressure become 

m = -ll— _ ^^)V^^ + (^2 - 0^)^22 + ^3V^33]. 



(52) 
(53) 

(with a obvious modification 
(54) 



Using this relation in Eq. ( p6| ) we find 

+ (7 + 1)^3] F33 + (7 - 1) [(1 - ^''Wii + {A2 - n')V22\ = 0. 

Equations (p6|), (p7|), ( p9| ) and (p2|),(|5^) can be manipulated to the the following set 

[A^ + 2 (1 - ^f)] Vii - [A^ + 2 {A2 - ^'^)] V22 - AXW12 = 0, 
[a2 + 2 (1 - n')] Vix + [A^ + 2 {A2 - n^)] V22 - 2(A2 + 2^3)1/33 + 4Ar! {Vi,2 - ¥2,1) = 0, 

(A^ + 1 + ^2 - 2^(2) 1/12 + ^MVii -¥22)= 0, 

X^{Vi,2 - 1^2;l) - nX{Vii + 1/22) + (1 - A2)Vi2 - 0, 



(55) 

(56) 
(57) 
(58) 
(59) 

which should be supplemented with Eq. (|5^). The corresponding characteristic equation is of order 8 and has been 
solved numerically. The results for the real parts of the modes (which are purely real) are shown in Fig. 4 as a 
function of A^ for several values of Q, for andx fixed ratio Vi^ /A2 = 0.1. 



V. CONCLUSIONS 

We have analyized the hydrodynamic oscillations of Bose-condensed atomic clouds at zero temperature in the 
Thomas-Fermi approximation. The equilibrium shape of the cloud in a axisymmetric trap represent either a prolate 
or oblate spheroid of revolution, which, for a particular choice of the rotation rate and trap potential, degenerates into 
a rotating sphere. The rotation frequency of the condensate is bounded from above by the characteristic frequency 
of the harmonic oscillator in a given trap (17^ /wq < 1). We have also analyized non-axisymmetric, triaxial ellipsoidal 
figures, which admit equilibrium solutions under additional constrain on the deformation in the plane orthogonal to 
the rotation axis A2 > Vi^ /loq. 

Small amplitude oscillations have been derived for linear perturbations from the rigidly rotating equilibrium back- 
ground state. The oscillations in axisymmetric traps, which are related to the transverse-shear and toroidal modes are 
found to be stable for all values of the trap deformation and its rotation frequency. These modes represent undamped 
oscillations (all eigenfrequencies are real) in the absence of dissipation, i.e., they are dynamically stable. A dynamical 
instability against the quasi-radial pulsation mode can arise when 7 < (1 — 31)^)/ (3 — $7^) < 1/3. Otherwise the 
system is stable, including the limit of the incompressible superfluid (7 — )■ 00), independent of the shape of equilibria. 
For a Bose gas with 7 = 2 the above stability condition is satisfied; hence, these modes represent stable oscillations 
independent of the rotation rate and deformation. Numerical results for the oscillations in non-axisymmetrical traps 
show that, both, the even and odd partity second harmonic modes of oscillations of triaxial clouds are stable for the 
adiabatic index 7 = 2 relevant to the interacting dilute Bose gas. 

Our results for axisymmetric traps (in which case analytical results are available) reduce to those of Ref. when 
the rotation frequency is set to zero and 7 = 2. Indeed, from Eqs. (22), (|^ and (|45|), which describe to the oscillations 
with ni = ±1, m = ±2 and m = 0, respectively, for I = 2, we find 



(t{1 = 2,m = ±1) = 

a{l = 2,m = ±2) = ±^/2, 



A. 



OA 1 

a(? = 2,m = 0) = 2 + :^±2 



16-16^3 + 9^2, 



(60) 
(61) 

(62) 
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which coincide with Eqs. (22), (23) and (24) of Ref. on making the appropriate changes in the parameters 
describing the trapping potential. 

It is instructive also to compare the results above to the classical analysis of the equilibrium and stability of the self- 
gravitating fluids |l9[| in view of speculations [^2| that intense off-resonant laser beams can give rise to a gravitational- 
type potential between the condensate particles leading to self-bound configurations. Self-gravitating systems, in 
particular the axisymmetric figures, are stable against transverse-shear and pulsation modes in the incompressible 
limit; the same is found in the above model of trapped rotating condensates. However, in the compressible case, the 
self- gravitating fluids are unstable against the pulsation modes whenever the adiabatic index 7 < 4/3 (the precise 
value of the critical 7 depends on the rotation rate) and are stable otherwise. In the present case, the system is 
unstable when 7 < (1 — — fi^), where the rotation frequency covers the range ^ci < < 1- Another major 

difference between the two systems is their stability against the toroidal modes. The axisymmetric self-gravitating 
fluids are unstable dynamically (i.e. in the absence of dissipation) beyond the point a§/aj = 0.05, as the deformation 
is increased. At a smaller deformation a\/a\ = 0.19, these oscillations become neutral, which is a prerequisite of the 
onset of secular (i.e. driven by the viscosity) instability. In contrast, the toroidal modes of the trapped condensates 
are always dynamically stable and there are no neutral points within the allowed parameter space where equilibrium 
figures exist. 

The present model can be extended for a study of the higher order {I > 2) harmonic oscillations by constructing 
higher order virial equations as well as to finite temperatures, in which case the viscosity of thermal excitation, hence 
secular instabilities, and mutual friction between the vortex lattice state and the excitations play a roleQ 
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APPENDIX: HYDRODYNAMIC EQUATIONS OF MOTION FOR ROTATING CONDENSATE 

As is well-known, dilute Bose gases can be described in the mean-field approximation in terms of the Gross-Pitaevskii 
(G-P) theory. The latter theory (in analogy to the phenomenological theories of superconductivity and of superfluid 
He^ near the A-point), is based on a Ginzburg-Landau type functional for the wave- function tp of the coherent state, 
whose variation provides the equation of motion for ij). In the stationary case this functional for a dilute Bose gas has 
the well-known form 



2m 



|V^(x)|^ 



1. 



-0tr(x)|V'(x)|^ + -{/o|V^(x)| 



cfx, 



(Al) 



where Uq = Anh^a/m. The G-P equation is obtained by taking functional derivative with respect to ip* subject to 
the constarint that the particle number N is constant. The extremum condition 6(E — fiN)/6'ijj* — gives 



^AV'(x) + \ [m(/.t,(x) C/o|^(x)n V(x) = /iV(x), 
2m 2 '■ ■' 



(A2) 



where the Lagrangian multiplier [i has the meaning of the chemical potential of particles. The time-dependent 
generalization of ( [A2| ) follows on the assumption that the temporal variations of il) should be described by a first order 
equation which, by analogy with the quantum mechanics, is written as 



^av'(x,i) ^ , 



y(/.t,(x)+C/o|^(x,t)pJ ^(x,t). 



(A3) 



'^Such a program for self-gravitating superfluids has been carried out recently in Ref. 
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On writing t) = r]{x, t)e^'^^^'*\ the superfluid density and velocity can be expressed as 



v(x,t) = —Vip(K,t). 
m 



(A4) 



The real part of Eq. (A3) after deviding by mrj, applying a gradient, and multiplying the result by p becomes 



Ar; 1 1 V , 

2m^ rj 2, m 



(A5) 



which is the Euler equation for the condensate in the zero-temperature, inviscid limit. It differs from the analogous 
equation for the ordinary fluids by the "quantum pressure" term oc A77/77. Note that any constant term can be added 
to the second bracket, for example, the chemical potential in the ground state. 

The imaginary part of Eq. (A3), on multiplying by 2mr]/h, leads to the mass conservation equation 



d_ 

di 



p + Vpv = 0. 



(A6) 



When gradients of the -0 are small, the "quantum pressure" term in Eq. ( A5) can be dropped and it reduces to 



d_ 

di 



vV 



-Vp- 



(A7) 



where p = {Uo/2m)p'^ = (2t: fi^ a /m'^)p^. Note that formally Eq. (A?) is identical to the Euler equation in the ordinary 
hydrodynamics; the distinctive feature of the superfluid is that the superflow is irrotational V x v = in general (the 
special case of the rotating superfluid, when the analogy becomes perfect, is discussed below). 
Eq. ( |A7| ) can be derived also starting from the momentum conservation equation: 



y- [0(x,t)VV*(x,t) - V*(x,i)V^(x,i)] 



where 



d 



(A8) 



Am? 



dip dip* 



dxi dxk dxidxk 



+ pSik, 



(A9) 



p is the pressure, c.c. stands for complex conjugate. The right-han d-sid e of E q. (A8) is the external force per unit 
volume. On writing ■(/'(x, t) = r]{x,t)e'''^^^'''^\ and using the relations (A4) Eq. (|A8|) becomes 



d_ 
di 



pvi 



Tp- [pViVk +pSik) = p^Vi^itr- 

dxk 2 



(AlO) 



The time derivative of p can be eliminated in terms of Eq. ( A6). If we use for the pressure of the dilute Bose gas the 
relation 



2TThap'^ 
m? 



64 f pa^ 



1/2- 



(All) 



and keep the leading order term in the diluteness parameter, we arrive again at Eq. ( A7). 

If the condensate is rotating at a angular velocity which is lar ger than the critical one f2ci , its energy is minimized 
via creation of vortices; then the curl of the last relation in (A_4) is non-zero, rather the phase of the superfluid order 
parameter changes by 27r around a path which encircles vortex lines 



m ' ' TO ^ 



(A12) 



where u is a unit vector along a vortex line, xj is the radius vector of a vortex line in the plane orthogonal to the 
vector v, and 5'-^^ is a two-dimensional delta function in this plane. If f2 ^ the macroscopic hydrodynamic 
equations involve only course-grained quantities, which are averages over a l arge number of vortices (i.e. over scales 
much larger than the size of a single vortex). The right-hand side of Eq. (A12) then becomes proportional to the 
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density of the vortex hnes, riy, since the continuum limit of vortex distribution implies (5^^-'(x — Xj) = riy. The 
left-hand side of Eq. (A12) in the course-grained limit gives since the energy is minimized by a superflow which 
mimics a rigid body rotation (this minimization is carried out e.g. , in e.g. Ref. p4|| ). Writing Eq. (A7) in the 
frame rotating uniformly with the angular velocity amounts to adding to the right-hand side of this equation the 
centrifugal potential \ fl x xp/2 and the Coriolis acceleration 2u x $7 (here u is the superfluid velocity in the rotating 
frame). With this substitution we recover Eq. (|l|). Note that the analogy to the Euler equation for a uniformly 
rotating ordinary fluid now is complete. 
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Fig 3: The real parts of the odd parity modes in non-axisymmetric traps for the fixed ratio AijVi?' = 0.1. Conventions are the same as in 
Fig. 1. 
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Fig 4; The real parts of the even parity modes in non-axisymmetric traps for the fixed ratio A2/f2^ =0.1 and 7 = 2. Conventions are the 
same as in Fig. 1. 
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